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Abstract 

This paper is concerned with the strong solution to the Cauchy-Dirichlet problem 
for backward stochastic partial differential equations of parabolic type. Existence 
and uniqueness theorems are obtained, due to an application of the continuation 
method under fairly weak conditions on variable coefficients and C 2 domains. The 
problem is also considered in weighted Sobolev spaces which allow the derivatives of 
the solutions to blow up near the boundary. As applications, a comparison theorem 
is obtained and the semi- linear equation is discussed in the C 2 domain. 
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1 Introduction 

m 

In this paper we consider the Cauchy-Dirichlet problem for backward stochastic partial 
CSj ■ different equations (BSPDEs, for short) either in the non- divergence form: 

dp(t,x) =— \a l3 (t,x)p x i x j(t, x) + b' l (t, x)p x i(t, x) — c(t, x)p(t, x) 

+ a lk {t, x)q*i(t, x) + u k {t, x)q k {t, x) + F(t, x)]dt (1.1) 
+ q k (t, x)dW k } (t, x) G [0, T) x V 



O 
O 



^ ■ or in the divergence form: 

h ; 

dp = - [{a l3 p x3 + a ik q k ) x i + b l p x , - cp + v k q k + F] dt + q k dW k (1.2) 
with the terminal-boundary condition: 

P (t,x) = o, te [o,T], xe&D, 

p(T, x) = (j)(x), x £ V. 



Supported by NSFC Grant #10325101, by Basic Research Program of China (973 Program) 
Grant # 2007CB814904, by the Science Foundation of the Ministry of Education of China Grant 
#200900071110001, and by WCU (World Class University) Program through the Korea Science and 
Engineering Foundation funded by the Ministry of Education, Science and Technology (R31-2009-000- 
20007). 

2 Department of Finance and Control Sciences, School of Mathematical Sciences, Fudan University, 
Shanghai 200433, China. E-mail: kduOfudan. edu. cn (Kai Du), sjtangOf udan. edu. cn (Shanjian Tang). 

3 Graduate Department of Financial Engineering, Ajou University, San 5, Woncheon-dong, Yeongtong- 
gu, Suwon, 443-749, Korea. 



1 



Here D is a domain of the ci-dimensional Euclidean space, and W = {W 7 ^; t > 0} is a 
G?i-dimensional standard Wiener process, whose natural augmented filtration is denoted 
by {^t}t>o- The coefficients a,b,c,a,u and the free term F and the terminal condition 
are all random fields. An adapted solution of equation (11. ip or (11. 2p is an 8? x B(T>)- 
measurable function pair (p,q), which satisfies, in addition to ( II. 3p . equations (11. ip or 
( 11.21) under some appropriate sense, where 2? is the predictable cx-algebra generated by 
{^t}t>o- 

BSPDEs, which are a mathematically natural extension of backward SDEs (see e.g. 
[3 E]), arise in many applications of probability theory and stochastic processes, for 
instance, in the optimal control of SDEs with incomplete information or more generally of 
stochastic parabolic PDEs, as adjoint equations (usually in the form of (11. 2p ) of Duncan- 
Mortensen-Zakai filtering equations (see e.g. (2j [161 EE EQl [25] ) to formulate the stochastic 
maximum principle for the optimal control, and in the formulation of the stochastic 
Feynman-Kac formula (see e.g. [TJ]) in mathematical finance. A class of fully nonlinear 
BSPDEs, the so-called backward stochastic Hamilton- Jacobi-Bellman equations, appears 
naturally in the dynamic programming theory of controlled non-Markovian processes (see 
HE]). For more aspects of BSPDEs, we refer to e.g. P [2Tj [221 123] • 

Equation ( II. 2p is usually understood in the weak sense (see Definition ^. II (ii) in Section 
2). When the coefficients a and a are differentiate in x, equation ( II. ip can be written into 
the divergence form (11. 2p . The existence and uniqueness of the weak solution of equation 
( 11. 2p in the whole space M. d follows from that of backward stochastic evolution equations 
in Hilbert spaces (see e.g. [5l Prop. 3.2]). However, weak solutions have low regularity, 
which find difficulty in many applications. Strong solutions and even classical solutions 
are required in many occasions. In the case of T> = R d , the theory of strong solutions on 
BSPDEs is now fairly complete. For instance, a H^-theory of the Cauchy problem for 
BSPDEs can be found in [21 EJ [111 EE 121] • On the contrary, there are very few discussions 
on the Cauchy- Dirichlet problem for BSPDEs. Here we could mention only two works. A 
special form of equation (II. 2p with Dirichlet conditions is studied in [23] by the method 
of semigroups, in the context that the coefficients are independent of (w,t). The other is 
our previous work [6], where the equations are solved in weighted Sobolev spaces. 

A main difficulty in strong solution of the Cauchy-Dirichlet problem for BSPDEs (and 
SPDEs) is to estimate the second order partial derivatives of the solution. In the theory of 
deterministic PDEs, it is solved with the help of the estimate of the derivative in t, which 
makes any sense in general neither for BSPDEs nor for SPDEs. For SPDEs, Flandoli pi)] 
establishes some regularity under additional compatibility conditions, and Krylov [12] 
studies the equations in weighted Sobolev spaces allowing the derivatives of the solutions 
to blow up near the boundary of T>. Note that there is an essential difference between 
SPDEs and BSPDEs: the noises in the former are exogenous and play an active role, while 
in the latter they are governed by the randomness of the coefficients and the terminal 
condition and thus they are endogenous, coming from martingale representations. The 
regularity of BSPDEs turns out to be more like deterministic PDEs than that of SPDEs. 

In this paper, we prove the existence and uniqueness of the strong solution (see Defini- 
tion [27TJ (i) in Section 2) of equation (II. ip without involving any additional compatibility 
condition nor any weighting functions. Our approach is based on the method of the 
odd reflection and some classical techniques from the theory of deterministic PDEs. Our 
assumptions on the coefficients are rather natural since they include the rather general 
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deterministic PDEs where the leading coefficients are not necessarily differentiable in 
x. Unfortunately, in contrast to deterministic parabolic PDEs (see e.g. Theorem 7.1.6 
in Evans |9J), further regularity for BSPDEs seems to be hopeless since the unknown 
random fields are not expected to be differentiable with respect to t as in the classical 
sense. However, we can consider the equations in weighted Sobolev spaces which allow 
the derivatives of the solutions to blow up near the boundary. Starting from the existence 
and uniqueness of the strong solution, we prove a slightly different version of our previous 
work |6j, and obtain the interior regularity for equation (II. ip . In the last part of our 
paper, we prove a comparison theorem for the strong solution of equation (II. ip . and we 
also discuss a class of semi-linear BSPDEs in C 2 domains. 

The rest of the paper is organized as follows. In Section 2, we introduce some notations 
and preliminary results. In Section 3, we state our main existence and uniqueness result 
in Theorem 13. 11 on the basis of which we study the equations in weighted Sobolev spaces. 
Section 4 is devoted to the proof of Theorem I3.1[ which is divided into two subsections. 
Finally, in Section 5, we prove a comparison theorem, and discuss semi-linear BSPDEs in 
C 2 domains. 

2 Preliminaries 

2.1 Notations 

Let (f2, JP, {^t}t>o, P) be a complete filtered probability space on which is defined a d\- 
dimensional standard Wiener process W = {W t ; t > 0} such that {^t}t>o is the natural 
filtration generated by W, augmented by all the P-null sets in JF. Fix a positive number 
T. Denote by & the cx-algebra of predictable sets on Q x [0,T) associated with {J^t}t>o- 

Let T> be a domain in M. d with boundary of class C 2 . 

For the sake of convenience, we denote 

DiU = u x i , D^u = u x i x j , Du = u x = (Am, . . . , D d u), D 2 u = (Ajw)dxd, 
and for any multi-index a = («i, . . . , ot d ) 

D a = D* 1 ■ ■ ■ D a d d , \a\ = «! + •• - + a d . 
For any two multi- indices a = («i, . . . , c^) and (3 = (/3 1; . . . , f3 d ), we define 

a + p = («i a d + p d ). 

We shall also use the summation convention. 

Now we introduce some function spaces. For any integer k > 0, we denote by C k (T>) 
(or C k (T>)) the set of functions having all derivatives up to order k continuous in T> (or 
V); by C^iV) (or Cj;(V)) the set of those functions in C k (V) (or C k {V)) whose partial 
derivatives up to order k are uniformly bounded in T> (or T>). 

For a given Banach space £>, denote by L 2 (Q x (0, T), B) the space of all B-valued 
predictable process X : Q x [0, T] — >• B such that E J Q T \\X(t)\\ 2 B dt < oo. 
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Let H m (V) be the Sobolev space W m > 2 (V) and H™{V) = W m ' 2 (P). Denote 

H°(D) = L 2 (fi x (0, T), L 2 (D)), 

H m (X?) = L 2 {Q x (0, T), H m {V)), m = -1,1,2,..., 

HJ(P) = L 2 (fi x (0, T), flJ(P)), 7i = 1, 2, 3, ... , 

n h 2 (£>) = L 2 (fi x (o, T), ^, # X (£>) n # 2 (£>)). 

Note that = H 1 (IR d ), Hp n M 2 (R d ) = U 2 (R d ). In addition, denote 

|| • 1 1 0,2? = || ' IU 2 (I?)) II ' ||m,2? = || ' ||_ff m (2?), 771 = — 1, 1, 2, 

Moreover, for a function u defined on Q x [0,T] x D, we denote 

\h\f m ,v = E t \\u(t,-)\\ 2 miV dt, m = -l,0,l,2,.... 
Jo 

The same notations will be used for vector-valued and matrix-valued functions, and in 
the case we denote |w| 2 = ^ |ii*| 2 and |w| 2 = \ u ^\ 2 j respectively. 
Furthermore, we define the two product spaces 

H 2,1 (V) = (Ml n H 2 (£>)) x H 1 (P; R dl ), 

U 2 ' 1 (V) = (Hj HI 2 (£>)) xHj(D;l 4 ), ( "' U 
both being equipped with the norm 

||(u,i;)|| W 2,i(2j) = ( HI u \\\\ v + HI v \\\l v ) 1/2 . 
It is clear that both U 2 ^{V) and H Q ' (£>) are Banach spaces. 



2.2 An Ito formula 

Let V and H be two separable Hilbert spaces such that V is densely embedded in H. We 
identify H with its dual space, and denote by V* the dual of V. We have V C H C V*. 
Denote by || ■ \\h the norms of H, by (■,■)# the scalar product in i/, and by (-, •) the 
duality product between V and V*. 

Consider three processes v, m, and v*, defined on Q x [0, T], taking values in V, H and 
V*, respectively. Let v(co,t) be measurable with respect to (u,t) and be ^-measurable 
with respect to u for a.e. t. For any r\ G V, the quantity (r],v*(co,t)) is ^-measurable 
in a; for a.e. t and is measurable with respect to (u,t). Assume that m(u,t) is strongly 
continuous in t and ^-measurable with respect to u for any t, and that it is a local 
martingale. Let (m) be the increasing process in the Doob-Meyer Decomposition of the 
sub-martingale ||77i||jj (see e.g. [TSJ Page 1240]). 

Proceeding identically to the proof of Theorem 3.2 in Krylov and Rozovskii [T3], we 
have the following result concerning Ito's formula, which is the backward version of [T3"l 
Theorem 3.2]. 
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Lemma 2.1. Let G L 2 (fi, JP T , H) . Suppose that for every i] <E V and almost every 
(oj,t) G fi x [0, T], it holds that 

(V, v(t)) H = (v, v)h + j (v, v*(s))ds + (r), m(T) - m(t)) H - 

Then there exist a set fi' C fi s.t. P(tt') = 1 and a function h(t) with values in H such 
that 

(a) h(t) is ^t-measurable for any t G [0,T] and strongly continuous with respect to t 
for any uj, and h(t) = v(t) (in the space H) a.e. (u,t) G fix [0,T], and h(T) = ip for 
any uj G fi'; 

(b) for any uj G fi' and any t G [0, T], 

\\h{t)\\ 2 H =y\\ 2 H + 2 J (v{s),v*{s))ds + 2 J (h(s),dm(s)) H -(m) T +(m) t . 

2.3 Notions of solutions to BSPDEs 

Throughout this paper, we assume that the given functions 

a = (a ij ) : Qx[0,T}xV ^ S d , a = (a ik ) : fi x [0, T] x V -> R dxdl , 
6 = (6*) : Ox[0,T]xD4l J , i/= (u k ) : fix[0,T]xD4 R dl , 

c, F: fix[fl,r]xD->l 

are all ^ x i?(P)-measurable (5 d is the set of real symmetric d x d matrices), and the 
function (j) : fi x V — > R is j^r x £? (Immeasurable. 

Let us now turn to the notions of solutions to equations (11. ip and (11.21) . Throughout 
this subsection it will be supposed that the coefficients of our equations, i.e., the functions 
a, b, c, a and u, are all bounded. 

Definition 2.1. A pair of random fields {{p{uj, t, x), q(uj, t, x))\ (oj, t,x) G fi x [0, T] x R d } 
is called 

(i) a strong solution of equation (II .ip with the terminal-boundary condition (11.31) if 
(p,q) G H 2 > l (V) and p G C([0, T), L 2 (V)) (a.s.) such that for each t G [0,T] and a.s. 
uj G fi, it holds that 

p(t, x) =(f>(x) + / [a 13 (s,x)Dijp(s,x) + b' l (s,x)Dip(s,x) — c(s,x)p(s,x) 

_ T ( 2 - 2 ) 

+ a ik (s, x)Diq k (s, x) + u k (s, x)q k (s, x) + F(s, x)] ds - q k (s, x)dW k 

for almost every x G R d ; 

(ii) a weak solution of equation (ll.2p with the terminal-boundary condition (11.31) . if 
(p,q) G HJ(D) x H°(£>;R dl ) such that for every 77 G Cq°(P) and almost every (w,t) G 
fix [0,T], it holds that 

p(t,-)r)dx= f <f)T)dx+ j j [(a ij Dip + cr^q^DjT] 
v Jv Jt Jv ^2 3) 

+ (UDip -cp + v k q k + F)rj] dxdt - [ [ q k r]dxdW k . 

Jt Jv 
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It follows from Lemma 12.11 that the first component of the weak solution of equa- 
tion (11.21) has a continuous version in L 2 {V), i.e., p £ C([0, T], L 2 (V)) (a.s.). For the 
strong solution of equation (II .ip with the terminal-boundary condition (11.31) . we have the 
following 

Proposition 2.2. Let (p,q) £ "H 2 ' 1 ^). The following statements are equivalent: 

(i) (p, q) is a strong solution of BSPDE (11.11) and fll.3|) .' 

(ii) for every r\ £ C^(T>) and a.e. (u>,t) £ Q x [0, T], it holds that 



/ p(t, -)r}dx = / (f)7]dx + / [a^DyP + UDip — cp 
Jv Jv Jt Jv 

+ a ik Diq k + v k q k + F] rjdxdt - f f q k 7]dxdW t 

Jt Jv 



(2.4) 



k 



Proof. It is clear that (i) =>■ (ii). Now we prove (ii) =>■ (i). It follows from Lemma [27T1 that 
p £ C([0, T], L 2 (T>)) (a.s.). Then from the time continuity of the (stochastic) integrals, 
we know that equation ( 12.4ft holds almost surely for every 77 £ C^(T>) and all t £ [0, T]. 
On the other hand, since (p, q) £ "H 2,1 (£>), both sides of equation ( 12.21) (as functions 
of x) belong to L 2 (T>) a.s. for every t. Since C^(T>) is dense in L 2 (V), we know that 
equation (12. 2\i holds in the space L 2 (V) for every t, which evidently implies (i). The proof 
is complete. □ 

Remark 2.1. Note that the space of test functions C^{T>) in Definition 12.11 and Propo- 
sition [272] can be replaced with Hq(V). 

If the coefficients a and a are differentiable in x, then equation ( 11. ip can be written 
into the divergence form (ll.2p , which allows us to define the weak solution to (II. ip . Then 
Proposition 12.21 indicates that a strong solution of (II. ip is a weak solution of (II. ip . 

The following lemma concerns the existence and uniqueness of the weak solution of 
equation ( II. 2p . which is borrowed from Proposition 3.2 in [5]. On the other hand, it can 
be proved by the duality method as in Zhou [21] and Lemma 12.11 

Lemma 2.3. Assume that kI + aa* < 2a < for some constant k > and 

that the functions b l ,c and v k are bounded by K. Suppose that F £ Mr l (T>) and £ 
L 2 (Q, J^t, L 2 (T>)) . Then equation ( 11.2j) with the terminal-boundary condition ( 11.3j) has a 
unique weak solution (p,q) in the space Hj(X>) xM°(D; R dl ) such thatp £ C([0, T], L 2 (D)) (a.s.) 
and 

ll|p||l!i> + lllfflllo l i>+£ SU P lb(^-)llo^<C(lll^lll 2 i^+^H0llo^) 5 (2-5) 

0<t<T 

where the constant C = C(K, k,T). 



3 Existence, Uniqueness, and Regularity 

In this section, we state our main results on the existence, uniqueness and regularity of 
the strong solution of equation (II. ip with the terminal-boundary condition (11.31) . 
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3.1 Existence and uniqueness of strong solutions to BSPDEs 

The weak solution of a deterministic parabolic PDE can be shown to belong to the space 
L 2 (0, T; H 2 (V)) if the free term belongs to L 2 ((0,T)xD) and the initial data lies in H%(V) 
(see e.g. Evans [9]). Flandoli [10] formulates a counterpart for a SPDEs with additional 
compatibility conditions on the free term. In the following, we obtain a counterpart for a 
BSPDE, which, in a remarkable way, does not require any compatibility condition like a 
SPDE. The higher regularity of the (strong) solution allows us to weaken the assumptions 
on the leading coefficients a and a so that they are not necessarily differentiable in x, 
where equation (II. ip is difficult to be written into the divergence form (II. 2p . 

Fix some constants K G (1, oo) and po, k G (0, 1). Denote 

B + = {x G R d : \x\ < 1, x 1 > 0}, B p (x) = {y G R d : \x - y\ < p}. 

Assumption 3.1. For every x G &D there exist a domain U C B$k Po (x) containing the 
ball -E>4 Po (x) and a one-to-one map $ : 2B + — > U PI D having the properties: 

x = $(0), D B Apo (x) n V, <$>{dB + n {x 1 = 0}) C dV; 

^\ 2 <\(D^\ 2 <^\ 2 , V(el d ; 

< K for any multi-index a s.t. |a| < 2, 

where is the Jacobi matrix of $. 

Note that in view of the Heine-Borel theorem, Assumption 13.11 is true if the domain 
T> is bounded and its boundary is C 2 . 

Assumption 3.2. The super-parabolicity condition holds: 

kI + oo* <2a< V (w, t, x) G n x [0, T] x £>. 

Assumption 3.3. There exists a function 7 : [0, 00) — > [0, 00) such that (i) 7 is continuous 
and increasing, (ii) 7(7") = if and only if r = 0, and (iii) for any (00, t) G Q x [0, T] and 
any x, y G V, 

\a(u,t,x) - a(w,t,y)\ < 7(|x - y\), \a(u,t,x) - a(u,t,y)\ < j(\x - y\). (3.1) 

We have the following existence and uniqueness theorem, whose proof will be be given 
in the next section. 

Theorem 3.1. Let Assumptions [X71 HLM and \3.c>\ be satisfied. Assume that the functions 
b\c, and v k are bounded by K. If F G B°(V) and <p G L 2 (fi, & T , H%(V)), BSPDE (TTT]) 
and (11.31) has a unique strong solution (p,q) G T-Lq (T>) such that 

P GC([0,T],L 2 (P))nL oo ([0,T], J fJ 1 (I?)) (a.s.). (3.2) 

Moreover, we have the following estimate 

III P \\\l,v + III q \\\lv +E sup \\p(t, -)\\l v < C{ HI F \\\l v +^11011?^), (3.3) 

0<t<T 

where the constant C only depends on K,p , k,T, and the function 7. 

Remark 3.1. Since all constants C in this paper are independent of d\, all our results in 
this paper may be extended to the more general equation (11.11) where the di-dimensional 
standard Wiener process is replaced with a cylindrical Wiener process. 
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3.2 Solution in weighted Sobolev spaces and regularity 

Unfortunately, we could not establish any higher regularity for BSPDEs to correspond to 
the theory of deterministic parabolic PDEs, as given by Evans [9j Theorem 7.1.6], since 
the unknown functions are not expected to be differentiable with respect to t as in the 
deterministic sense. However, we can consider the equations in weighted Sobolev spaces 
allowing the derivatives of the solutions to blow up near the boundary, and furthermore 
obtain an interior regularity for BSPDE (11.11) and (11.31) . 

Let ip £ Cl(T>) be a nonnegative function such that ip(x) = for any x £ dT>. Fix a 
positive integer n. 

Assumption 3.4. For any multi-index a such that \a\ < n, we have 

i/j\<*\(\D a a ij \ + \D a b l \ + \D a c\ + \D a a ik \ + \D a v k \) < K, 
ijj\ a \D a F £ H°(£>), ^D a (j) £ L 2 (tt, & T , Hl{V)). 

Note that Assumption 13.41 implies the boundedness of the functions b, c and v, but 
does not imply Assumption 13.31 We have the following 

Theorem 3.2. Let Assumptions {3J\ El GDJ and\3J\ be satisfied. Then BSPDE (O]) 
and (11.3P has a unique strong solution (p,q) such that for any multi-index a s.t. \a\ < n, 

{^\Dy,^D«q) £%f{V), 

^H jD>GC '([0,T],L 2 (P))nL oo ([0,T],iJ 1 (P)) (a.s.), 

and moreover 

HI ijj M D a p \\\ 2 2V + HI ^D a q \\\ 2 lv +E sup \\^D a p{t,-)\\ 2 lv 

0<t<T 

<CJ2 (ill W\lv +m m D^\\l v ) , (3.5) 

l/9[<H ' 

where the constant C only depends on the norm of ip in C 2 (T>), the parameters K,po,n, 
and T , and the function 7. 

We need the following lemma, which can be found in [T2] . 

Lemma 3.3. // both v and ipDv lie in L 2 {V), then ifjv £ Hq(V). 

Proof. Define K n = {x £ T> : dist(x, &D) > 4 • 2~ n }. Take a nonnegative function ( £ 
Cg°(R d ) such that suppC C #i(0), J Rd ( = 1. Define ( n (x) = 2 nd ((2 n x) and r? n = 
C n * l Kn . We have supp(r? n ) C K. n+ i and r] n \ Kn _ 1 = 1. Since ip £ C 2 (V), we have \i>(x)\ < 
Cdist(x, dV). It is not hard to show that \r] n \ < 1, \ipDrj n \ < C, and r) n ipv £ Hq(V). Then 
we can get that both rj n ipv — > ipv and r] n D{%l)v) —> D(ipv) strongly in L 2 (V) as n — )■ 00, 
and moreover 

%^)-D(H| 2 <2 / |^ n |>| 2 + 2 / \ Vn D(iJv) - D(i;v)\ 2 
v Jv Jv 

<C I \v\ 2 + 2 [ \rj n D{ipv) — D(ipv)\ 2 

JV\K,„-i JV 

—7-0 as n —7- 00. 

Hence ipv £ H%(V). □ 
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Proof of Theorem \3.2i The proof consists of two steps. We suppose for the moment that 
ip G C n+2 (T>), which will finally be dispensed with. 

Step 1. We first prove that Theorem 13.21 is true for the case where the leading coef- 
ficients a and a are differentiable in the space variable x with the gradients a x and o x 
being bounded and thus equation (11.11) can be written into the divergence form (jl.2p . 

We use the induction. Theorem 13.11 shows that Theorem 13.21 is true for the case of 
n — 0. Assume that it is true for the case of n — m — 1 (m > 1). It is sufficient for us to 
show that it is true for the case of n = m. 

We assert that (ip m D a p,iJj m D a q) G x B°(V;R dl ) for any multi-index a s.t. 

\a\ = m. Indeed, we know from our assumption that for any multi-index /3 s.t. \/3\ < m—1, 

^m-i D i3 p G jja^p), iff^Dfq E H 1 (P; R dl ). 

Keeping in mind that ip G C%(V), we can easily get by induction that 

ij m - 1 D^p x G U^V), ^) m - x D^p xx G H°(D), and ^ m - l D^q x G H°(D; (3.6) 

Then we have 

^m-i D a p e c jjO^ ^D(i) m ~ 1 D a p) G H°(X>), and V m £> a <? G H°(£>; M ' 1 ). 

In view of Lemma [3.31 we have from the first two relations that ip m D a p G Hj(X>). 

Take any r? G C °°(P). Since ^ G C 6 n+2 (£>), we know that D a (ip m 7]) G #o(£>). From 
Proposition I2.2[ we have 



p(t,-)D a (^ m ri)dx- / <PD a (i) m ri)dx 



p 



v 



t Jv 

T 



a ij D ijP + UDip -cp + a ik D iq k + v k q k + F D a (ij m r])dxdt 

[ [ q k D a (i; m ri)dxdW k , a.e. (u,t) G tt x [0,T]. 

Using the integration by parts formula, we show that the function pair (ip m D a p, ip m D a q) G 
Mo(X>) x B°(V; R dl ) is a weak solution of the following BSPDE 



k 

t > 



du = - [a ij D ijU + a lk D lV k + F]dt + v k dW i 
u{t,x) =0, t G [0,T], x G <9£>, 
u(T,x) = {ip m D a (j)){x), xeV 

with u and t> being the unknown functions. Here 
F = [O 1 ^^') ('••" iyp, , ) + {^D^a lk ) {iJ hl D^q k t ) 

/3+-y=a,\/3\>l 

+ [{^ m D^)(^D^ Pxl ) - (^D?c){^D~<p) 

/3+7=a 



(3.7) 



(3. 



- 2ma ij ip m - 1 D i ^D a p x] - ma ij ip m - l D ij ^D a p - 



'j ' 

/3+7=a 



m(m 



l)a ij ^ m - 2 D i ^D j ^D a p 



From (I3.6P and Assumption I3.4[ we see that F G M°(V). Moreover, from estimate (I3.5P 
for n = m — 1 (as a consequence of the induction assumption), we have 



l-^lllo.D ^ C 



< C 



E ( HI * WDf, P W\h> + HI HI?,* ) + \U m D a F lllo. 



V 



\<m-l 



£ ( HI ^ D ^ F MSU +E\\^D^\\l^ + \U m D a F 

\P\<m-\ 



|2 

l0,£> 



where the constant C only depends on the norm of tp in C 2 (V), the constants K,p ,K 
and T, and the function 7. Note that if) m D a <f) G Hq(P). Therefore, applying Theorem 
13. II to BSPDE (13.71) . we have Theorem 13.21 for n = m. 

Step 2. Now we remove the boundedness assumption on a x and a x made in Step 1. 

In view of Theorem 13.11 BSPDE (11.11) and (II. 3p has a unique strong solution (p,q). 
Due to Assumption 13.31 we can construct (e.g., by the standard technique of mollification) 
two sequences a r and a r with their first-order derivatives in x being bounded, which 
converge uniformly (w.r.t. (00, t, x)) to a and a, respectively, such that all a r and a r 
satisfy assumptions 13. 2\ I3.3[ and I3.4| with k in assumption 13.31 being replaced with k 2 . 
Then from Theorem I3.1| the following equation (for each n) 



dp r = - 

Pr\x£dV 



(ai j D ijPr + bW^ - cp r + afD { q k r + v k q k r +F)dt + q k dW i 



t j 



Pr\t=T 



has a unique strong solution (p r ,q r ) G 1-t^ l (V), which satisfies estimate (I3.3P with the 
constant C being independent of r. Then we can check that {(p r , q r )} is a Cauchy sequence 
in the space T-Lq (T>), whose limit is (p, q). Similarly, we have that {(tfj^D a p r , tjj^D a q r )} 
(| a | < n) is also a Cauchy sequence in the space 1-L^ l (T>), whose limit is denoted by 
(u a ,v a ). Evidently, we have that u a G C([0, T], L 2 (V)) n L°°([0, T], tf 1 (£>)), and 



[U a , V a ) ||^2 



1{V) + E sup \\u a (t,.)\\l v <C T ( \H W D^F\\\l v +E\\^D^\\l v ). 

(3.9) 

On the other hand, for every rj G Cg°(X>) and a.e. {oo,t), we have (|a| < n) 

(i)\ a \D a p r) r]) = (-l) |Q| (j9 r ,D Q (^ |Q| 77)) ->• (-l) |Q| (p,D Q (^ |a| r/)) = {^ a \D a p,rf), 

where we denote (u,v) = J v u(x)v(x)dx. Thus ip^D a p = u a . Similarly, we have 
ijr a]{ D a q = v a . Estimate (13.51) is derived from inequality (13. 9p . 

It remains to remove the additional assumption made at the beginning that ip G 
C^ +2 (V). Note that the constant C in our estimate only depends on \ip\c^ (and other 
parameters), which allows us to approximate tp in C 2 (V) by a sequence of nonnegative 
C^ +2 (P)-functions vanishing on the boundary. Then in view of Lebesgue's dominated 
convergence theorem, the proof is complete. □ 

By choosing a proper weighting function ip, we obtain the following interior spacial 
regularity for the strong solution of BSPDE (II. ip and (11.31) . 
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Corollary 3.4. Let Assumptions \3.1 \ and \3.& be satisfied. In addition, suppose that 
{\D a a ij \ + |D a 6*| + |D a c| + \D a a ik \ + \D a u k \) < K, 

\a\<n (3.10) 

FGH"(D), G L 2 {Q, & t , Hl{V) n H n+1 (V)). 

Here the integer n > 1. Tnen BSPDE (11. ip and (II. 3p nas a unique strong solution 
(p,q) G "Hq' 1 ^) snc/i tnat 

(i) the functions p and q satisfy (I3.6P and (13.31) ; 

(ii) p G H^+ 2 (P), g G e^DjR*), p G C([0, T], H™ oc {V)){a.s.), i.e., for any domain 
V CC V, we have 



V \\\l+2,V> + III Q \\\l+l,V> + E SU P IIp(*> 0IIS+1.2J' 

0<t<T 



(3.11) 



< C £ (p A l)- 2 (— )( HI F |||^ +F||0||^ +li23 ), 

m=— 1 

wi/i p = dist(X>', &D) and with the constant C = C(K, po, k, T); 
(iii) moreover, if n — d/2 > 2, then 

peL 2 (ttx (0, T), 0>, C\V')) n L 2 (fi, C([0, T) x £>'))> 
g G L\Q x (CT^^C^P'jR*)), 

/or any domain V CC X>. 

Proof. In view of Theorem 13 .![ it remains to prove the assertions (ii) and (iii). 
Without loss of generality, we suppose p < 1. Define 

K = {x G V : dist(x, <9D) > p/2}. 

It is clear that T>' C /C. Take a nonnegative function £ G Cg°(R d ) such that 

sup P Cci?£(0), / C = i, I^CI^Cp- 1 , l^d^Cp- 2 . 

Define ip = (p/2)£ * 1/e- It is not hard to show that ip is a well defined weight function 



for Theorem 13.21 and moreover 

In view of (13. 8p and keeping in mind (I3.12p . we show that for n = 1, 



^P|||L?+lll^lllL>+£ SU P llV^frOII? 



V 

0<t<T 



< C(p 2 HI D x p \\\l v + HI p ||| 2 iP + HI g \\\l v + HI ^DF Hl 2 ^ + HI i)D<\> \f i;D ), 
and for 2 < m < n and multi-index a s.t. |a| = m, 

l D a p\\\l v + \\\^ m D a q\\\l v +E sup U m D a p(t,-)\\l v 



0<t<T 

< c J2 ( HI ^ m ° p p W\h> + 111^^111^) 



|<ro-l 



m ^F||| 2 +F||^ m J D>|| 2 
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By induction, we have (here \a\ = n) 

\H n D a p\\\l v + \\\rD a q\\\l v +E snp U n D a p(t, ■) \\{ v 

0<t<T 

< Cp-'lWDxplWlv+cY; ( III III ^0111^)- 

|j8[<n 

By multiplying p~ 2n on both sides, we can easily get (13.111) . The assertion (iii) follows 
from the Sobolev embedding theorem. The proof is complete. □ 

Remark 3.2. In the case of n — d/2 > 2, the function pair (p, q) satisfies equation ( 12. 21) 
for all (t, x) G [0, T]xD and wgQ' s.t. F(Q') = 1, which is a classical solution of equation 
(II. ip with the terminal-boundary condition ( II. 3p (see e.g. [T5] for details). 



4 Proof of Theorem 3.1 



The proof is rather long, and it is divided into two subsections. We first consider the 
simpler domain of a half space, and then go to the general C 2 domain. 

4.1 The case of the half space 

In this subsection, we consider BSPDE ( II. ip and (ll.3p living in a half space. 
Recall R d + = {x G R d : x x > 0}. Denote y = (x 2 , ...,x d ). 

Definition 4.1. We say a function / defined on M. d has the property of reflection invari- 
ance, if f(x 1 ,y) = — /(— x x ,y) for almost every (x x ,y) G W 1 . 

For a function u defined on R^, define u and u as follows: 



u, on {x 1 > 0}, _ ( i 

0, on {x 1 < 0}; 



u{x l ,y), if x 1 > 0, 

u{x\y)= { 0, if ^ = 0, (4.1) 

— u(— x x ,y), if x 1 < 0. 



It is clear that u has the property of reflection invariance. 

Lemma 4.1. (a) Let m be a positive number. Then a function u G /f™(lR^) if and only 
ifu G H m (R d ). 

(b) The function u G i^lM */ and onl V % f^E E X (^ A ). 

Proof. The proof of assertion (a) can be found in, e.g., Chen [3j Page 48]. The necessity 
of assertion (b) follows from assertion (a). It remains to prove the sufficiency of (b). 
Indeed, we can find ip n G C^°(lR a! ) such that ip n — )■ u in if 1 (IR d ) as n — > oo. Denote 
Tp^x 1 ^) = — (p n (— x 1 , y). Note that u(x x ,y) = — u(— x x ,y). Thus we have Tp n — y u in 
J ff 1 (M d ). Define i\) n = [ip n + Tp n )/2. Then ip n -> u in H x {M. d ). Since ip n (0,y) = 0, the 
restriction of ip n on Mi belongs to if^M^). Thus u G H 1 (Mi). The proof is complete. □ 

The following existence and uniqueness result concerning the Cauchy problem of 
BSPDEs is taken from Du and Meng [51 Prop. 4.1], which can also be proved by means 
of the duality method of Zhou |24J and Lemma 12.11 
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Lemma 4.2. Consider the following BSPDE (on M. d ) 

dp = -[a ij (t)D ijP + a ik (t)D iq k + F]dt + q k dW t k 



p(T, x) = <fi(x), x G 



(4-2) 



where a and a are two predictable processes taken values in S n and M. dxdl , respectively, 
such that kI + aa* < 2a < kT 1 /, V(w,t). Suppose F G H°(R d ), G L 2 (ft, J? T> tf 1 ^)). 
T/ien BSPDE (121) fcos a urague strong solution (p,q) in M 2 (R d ) x H 1 (M rf ;R dl ) such that 
p G C([0, T], i/ 1 (IR d )) (a.s.) ; and moreover, 

III P 111^ + III 9 III?,** +E SU P UK*, OIIU" ^ C '(«» T )( HI F W\U* + E W<f>\\l**)- ( 4 - 3 ) 

0<t<T 

Now we have the following 

Lemma 4.3. Let Assumptions \3.& be satisfied with D = R+. Assume that a and a are 
invariant in the space variable x. Suppose that F G H^R 1 ^) and G L 2 (Q, Hq 0^+))- 
Then the following BSPDE 

dp = -{a^^Dijp + a ik (t)D iq k + F]dt + a fe dW t fe , 

= 0, x G <9M^, (4.4) 
p(7» = 0(x), I6lf 

has a unique strong solution (p,q) G %o' 1 (R+) swc/i that p G C([0, T], if 1 (R+)) (a.s.), and 
\\(p,q)\\ 2 H2 , 1{Ri) +E sup b(t,-)||? Rl <C(|||F|||2 Ml +E||0||2 R A (4.5) 

0<*<T + + + 

where the constant C depends only on k and T. 

Proof. Recalling the definition fl4TT]) . we have F G e°(R d ) and G L 2 (ft, & T , H x (R d )). 
From Lemma i2l the BSPDE 

dP = -[aVtyDyP + or ifc (t)AQ fc + F)dt + Q k dW k , . . 

P(l» = 0(x), a; G R d 1 ' ' 

has a unique strong solution (P,Q) in the space H 2 (R d ) x EI 1 (R d ;R dl ) such that P G 
C([0,T],H 1 (R d )) (a.s.), with the estimate 

III P + III Q III?,- sup ||P(t, OH^d < C(k,T)( HI F ||| 2 iRd +F||0|| 2 jRd ). (4.7) 

0<t<T 

By symmetry and the uniqueness of the solution (of equation (14.61) ). we know that P and Q 
have the property of reflection invariance, for a.e. (u, t). Denote by p and q the restrictions 
of P and Q on R+, respectively. From Lemma [4.11 (b), we know that p G Hq n H 2 (R c f_) 
and q G Hj(R^;R dl ). Moreover, p G C([0, T], If X (R+)) (a.s.). It is evident that the pair 
(p, q) is a strong solution of equation (14.4p . Since every strong solution of equation (14.4p 
is also a weak solution, from the uniqueness of the weak solution, we know that (p, q) 
is the unique strong solution of (14.41) . Estimate (I4.5p follows from inequality (14.71) . The 
proof is complete. □ 
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Now we prove the following perturbation result, which will be used in the proof of 
Theorem 13.11 

Proposition 4.4. Consider the following BSPDE 

dp = -[a ij D ijP + a ik D iq k + F]dt + q k dW t k , 

p(t,x) = 0, x E dR d +1 (4.8) 
p(T,x) = (j)(x), x E R d + . 

Assume that for a constant 5 > and for any (u, t, x) we have 

\a(t, x) - a (t) | < S, \a(t, x) - a (t) | < 5, (4.9) 

where {ao(t) : < t < T} and {cr (t) : < t < T} are predicable processes satisfying 
Assumptions^^, Suppose that F E B°(R d + ) and E L 2 (0, & T , H^(R d + )). 

Under the above assumptions, we assert that there exists a constant S(k,T) > such 
that if 5 < 5(k,T) then BSPDE (14.81) has a unique strong solution (p,q) E 'H 2 ) ,1 (R d h ) such 
that p E C([0, T],H 1 (R d _)) (a.s.), and 

\\(p,q)\\ 2 n2 , 1{K) + E sup |b(tr)lll Rl <^(lll^llloM+^llCn)' ( 41 °) 

0<t<T + + + 

where the constant C depends on k and T. 

Proof. In view of Lemma [4.31 we know that for any (u,v) E "H 2,1 (M^.), the BSPDE 

dp = — [do E>ijP + cT l Diq + (a y — cl^D^u 

+ (a i - 4)D iV + F]dt + qdW t , , . 

p(t,x) = 0, x E dR d + , 1 ' ' 

p(T,x) = <f>(x), xER d + 

has a unique strong solution (p, q) E Ul' l (R d + ) such that p E C([0, T],H 1 (R d _)) (a.s.). We 
define the operator A : "H 2,1 (R^) ->• n 2 '\R d + ) as follows: 

A(u,v) = (p,q). 

Then from estimate (14.51) . we obtain that for (ui,Vi) E 'H 2 ' 1 (R d _), i = 1,2, 

\\A{U! -U 2 ,V 1 -V 2 )\\ 2 H 2,i {R d +) < CSWim ~U 2 ,V 1 - V 2 )\\l i2 , 1{R dy 



Taking 5 = (2C)- 1 = (2C(k, T)) _1 , we have that the operator A is a contraction in 
% 2,1 (IR+), and thus there exists a unique element (p,q) E "H 2,1 (R+) such that A(p, q) = 
(p,q). Furthermore, we have (p, q) E ^(R^) and V e C([0, T], ^(Rj.)) (a.s.). It is 
clear that (p, q) is the unique strong solution of BSPDE (14. 8p . 

To establish estimate (I4.10p . in view of Lemma |4"73"| applying estimate (14.51) to equation 
dUTT]), we have 



^2,i(Mi) + E sup \\p(t, -)\\l Rd < C5\\(p,q)\\ 2 H2A{Ri) + C( HI F \\\l K +E\\(p\\l K ). 



0<t<T ' + y +' ' + ' + 

Taking 5 = (2C)" 1 , we obtain (14. lUp . The proof is complete. □ 



14 



4.2 The case of the general C 2 domain 

In this subsection, we shall complete the proof of Theorem 13.11 To simplify the notation, 
we define 

&\V) = {(u,v)en 2 \V):ueL\Q,L^([0,T},H\V))), 

ueC([0,T},L\V)) (a.s.)} 

being equipped with the norm 

H(^)iU^) = (\\(u,v)\\ 2 H2A{T)) + E sup \Ht,-)\\l v ) 1/2 . (4.12) 

0<t<T 

It is clear that S)^ l (V) is a Banach space. 
First we have the following fact. 

Lemma 4.5. Let $ be the map defined in Assumption \3.1\ and ^ be the inverse of $. 
Suppose u G #o( R +) s - 1 - supp(w) C B + U dR d + . Then mo$6 H%(D). 

Proof. The proof is straightforward. Take u n G C^°(B + ) such that u n — > u strongly in 
iJ 1 (Ml) as n — > oo. From the properties of it is easy to show that u n o $ g C 2 (T>) and 
supp(u n o$)cf/UP and then u n o \1> g Hq(V), where U is taken from Assumption 13.11 
Now we have 

\\u n o ^ - u o *||g )2J < | det(D$)| L o || Mn - 0, 

\\D(u n o ^ - u o < | det( J D$)| L o | J D^|| 00 || J DK - u)f Q)B+ 0, 

as n — ?■ oo, which implies mo$ g Hq(T>). □ 
The following is Theorem 13 . 1 1 under an additional assumption on the coefficients a and 

a. 

Proposition 4.6. Let the conditions of Theorem lff.il be satisfied. In addition, assume 
that a x and o~ x are bounded. Then BSPDE (II. ip and ( II. 3p has a unique strong solution 
{PA) e rf^iV) such that 

IIMU 2 2 ,2> + IIMIIi^ +£ sup \\p(t,-)\\l v <C{\\\F\\\l v +E\\m v ), (4.13) 

0<t<T 

where the constant C only depends K,p , k,T, and the function 7. 

Proof. Since a x and a x are bounded, equation (11.11) can be written into the divergence 
form 

dp = - [(a^p xJ + a ik q k ) x i + (V - a%)p xi - cp + {v k - a ik )q k + F]dt + q k dW k . 

From Lemma [2.31 BSPDE (11.11) and (11.31) has a unique weak solution (p,q) G Hj(X>) x 
H°(£>;M dl ). 

Now take a sufficiently small p G (0, po) to satisfy the following two conditions. 
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(1) j(8p) < 5 with the constant 5 = 5(k,T) being given by Proposition 14.41 In view of 
Assumption 13.31 for any (ou, t) and x, y G V, we have 

\a(t,x)-a{t,y)\<5, \a(t, x) - a(t, y)\ < 5 (4.14) 

if \x — y\ < 8p. 



(2) If x,y belong to the same domain U in Assumption 13.11 then for any (cu,t), 

\DV(x)a(t,x)(DV(x))* - DV(y)a(t,y)(DV(y)y \ < 5 U 

\D^(x)a{t,x) - DV(y)a(t,y) \ < 5x (4.15) 

if \x — y\ < 8p, where the constant 5i = S(k 2 , T) and \l/ is the inverse of $. 

Then take a nonneg ative function C G C^°(M d ) such that supp(C) C B 2p (0), ((x) = 1 
for \x\ < p. For any z G M. d , define for x G M. d , 

C(x) = ({x-z), p z {t,x)=p{t,x)C{x), q z {t,x)=q{t,x)C{x). (4.16) 
Then (p z , q z ) satisfies the following equation (in the sense of Definition 12.11 (ii)) 
dp* = - + a^qf + ( Z F + {b% z - 2a%* K< 

- (<>' J C,, + c)Cp + {v k - « lk C xi )Cq k ]dt + q z ' k dW k . 
In addition, define 

v z (x) = a x -^)- 

Now consider the following two cases. 

Case 1. dist(z, dV) < 2p . Then £>nsupp(p z ) C UfXD, where U is given in Assumption 
O Set x = $(y). Define 

u z (t, y) = p z (t, x), v z (t, y) = q z (t, x), (t, y) G (0, T) x R d y>+ . 

Obviously (u z ,v z ) G H^(M^ + ) x H (M^ + ;M dl ). Direct calculus shows that 

p z x r x s(t,x) = ¥ x r(x)¥ xs (x)u z viv j(t,y) + (Cp x i +Qip)(t,x)¥ x r x s(xWAy) 



(4.17) 



q x r k (t,x) = % r (x)v*f(t,y). 



Substituting the above relations into equation (I4.17p . it is not hard to check that the 
functions u z ,v z satisfy the BSPDE (in the sense of Definition 12.11 (ii)) 



du z = {n ij <i; nr + a ik v z /F)dt + v z > k dW k , 
w 2 |r^ + = 0, u z \ t=T = C 2 0, 

where (observe that u z = 0, v z = whenever r] z ^ 1) 
x = $(</), x = $(0), L° = a rs d 2 xrx s, 

a ij (t,y) = a rs (t,x)¥ xr %s(x)T] z (x) + a rs (t,x )¥ xr ¥ x s(x )(l -r, z {x)), 

a lk (t,y) = a rk (t,x)^)v z (x) + a rk (t,x )¥ xr (x )(l- V z (x)), 

F(t,y) = (CF)(t,x)+p x r(t,x)Q r 1 (t,y)+p(t,x)Q 2 (t,y) + q k (t,x)Q k (t,y), 

Ql(t,y) = (( z L°V)(t,xW yt (t,y) + (b r C-2a rs C)(t,x), 

Q 2 (t,y) = ((*rL°¥)(t,x)<S> r yi (t,y) - (L°( z + c( z )(t, x), 

Q k (t,y) = (u k C-a rk C z xr )(t,x). 



(4.18) 
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In order to apply Proposition 14.41 to BSPDE (|4.18p . we take 

a (t) = a(t, 0), a (t) = a(t, 0). (4.19) 

Note that supp(r^) C B ip (z). Then it follows from (14. 15|) that for any y G M. d and 
x = $(?/), we have (recall xo = $(0)) 

\a(t,y)-a{t,0)\ = \a rs (t, x)%M s (x) - a rs (t,x )% r ^Ux )\ ■ \ V z {x)\ 

< \a rs (t,x)% r ^is(x) - a rs (t,x )% r ^is(x ) \ ■ l B4p(z) 

< 5 1 = 5(k 2 ,T). 

Therefore, from Proposition 14.41 BSPDE (14. 18j) has a unique strong solution (u, v) such 
that 

u G Hq D EI 2 (IR^ + ), v G Hq(IR^ + ; M. dl ). 

It is clear that (u,v) is also a weak solution to BSPDE (I4.18p . From the uniqueness of 
the weak solution, we have u = u z and v = v z . Hence we deduce that 

u z G n H 2 (Mj i+ ), v z G R dl ). (4.20) 

Denote D(z,r) = B r (z) fl T>. Then (14.201) implies that restricted on the domain D(z,p), 
the solution 

(p, q) G U 2 (D(z, p)) x M L (D(z, p); R dl ) (4.21) 

for any z s.t. dist(z, &D) < 2p . 

Now applying estimate fl4TT0l to BSPDE (OHjl . we obtain that 

lib 2 + III IIIIr* + e su p K(v)ii? )R <i 

' «'+ ' «•+ 0<t<T y '+ 

< C(«,T)( III F |||^ + +E\\(C<t>) o $II? K+ )- 
On the other hand, it is evident that (recall supp(£ 2 ) C B2 P (z)) 

\\\P\\\2,D(*,P) + I" q HlW.P) +E SU P lb(*>0lll,D(*.P) 

0<t<T 

< \\\P Z \\\lv + \\\q z +Esup\\p z (t,.)\\l v 

t<T 

<C{\\\u z \\\l Rd +\\\v z \\\l d +E sup \\u z (t,-)\\l d ), 

V,+ v,+ 0<t<T v ' + 

II -^IHo.R 1 * , — ^( HI E \%,D{z,2p) + III P \\\l,D(z,2p) + III 1 \\\o,D(z,2p) )' 



^H(C^) o m K+ < CE\\C(f>\\l v < ceu\\I d(zM , 



where the constant C depends only on K and k. Therefore, we obtain 

III P \\\l,D{z,p) + III Q \\\i,D(z,p) + E q SUP^ \\p(t, 
< C( III ^ lllo,£>(z,2p) +^ll0lll,D(z,2/9) + lib llll,D(2,2p) + III Q \\\o,D(z,2p) )• 



(4.22) 



Case dist(z,<9X?) > 2p - This case can easily be reduced to the first one. Indeed, 
we can replace the domain T> by any half space with the boundary lying at a distance 
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2p from z. In this situation it is not necessary to flatten the boundary and to change 
coordinates. Then as above we deduce property (I4.2ip for any z £ V and obtain an 
estimate similar to ( I4.22p . 

Integrating both sides of inequality ( 14.221) over z £ M. d , we obtain that 

lib \\\lv + Ilk III + E su p lb(*r)lli,p 

o<t<T (423) 
< C( HI F 11^ +F||0||? )I? + ||b|||^+|||g|||^), 

where the constant C depends on K, p , k, T, and the function 7. Since (p, q) £ HIq(D) x 
H°(P;]R dl ), the right-hand side is finite. Recalling that (14.21 j) holds true for any z £ T>, 
the above estimate implies that the unique weak solution of BSPDE (II. ip and (11. 3p 
found by Lemma [2.31 belongs to the space "H 2,1 (D), and moreover, p £ C([0, T], L 2 (T>)) n 
L°°([0, T], (a.s.). From Proposition I2.2[ we know that (p, q) is the unique strong 

solution of BSPDE JTTJ and (PP . 

Replace the initial time zero by any s £ [0, T). Proceeding identically as before, we 
can obtain the following estimate similar to (I4.23P 



e[ \\p(t,-)\\l v dt + E [ \\q(t,-)\\l v dt + E sup |b(t, -)llf,x» 
< C (\\\F \\\l v +EU\\l v + E \\p(t, .) + jT ||g(t, Ollio* 



(4.24) 



In view of the definition of the strong solution (Definition 12. ip . we know that the pro- 
cess p{t, ■) is an L 2 (T>)-valued semi- martingale. Then applying Ito's formula for Hilbert- 
valued semi- martingales (see e.g. pQ, Page 105]), we have 

Ms, -)\\lv =H\\lv + 2 / T / piaVp** + Vp x > -cp + a lk q k t + v k q k + F)dxdt 

Js Jv 

\\q(t,-)\\l v dt-2 [ [ pq k dxdW t k . 

Js JV 



Taking expectations and using the Cauchy-Schwarz inequality, we have for any e > 

E [ T \\q{t,-)\\l,v < EU\\l v + 2E f [ p(a ij p xixj + b i p x , - cp + a ik q k xl + v k q k + F)dxdt 

Js Js Jt> 

< EU\\l v + eE I* {\\p(t,-)\\l,v+h(t,-)\\lv)dt 

J s 

+ C(e,K)E [ T \\p(t,-)\\lvdt+\\\F\\\l v . 

J s 

Taking e small enough and recalling (I4.24p . we have 

e[ \\p{t,-)\\l v dt + E I \\q{t,-)\\l v dt + E sup \\p{t,-)\\l v 

Js Js s<t<T 

< C( III F \\\l v +EU\\l v + E f \\p(t, -)\\lvdt), 

J s 



(4.25) 



where the constant C depends on K,p 0j k,T, and the function 7. In particular, we have 

E\\p(s, <c(\\\F +EU\\l v + e£ \\p(t, .)\\l v d?j . 
Using the Gronwall inequality, we have 

IIMIIi,* = [ T E\\p(s, -)\\lv < Ce CT { HI F +EU\\l v ). 
Jo 

The last inequality along with (I4.25P yields estimate (14.131) . 

It remains to prove q G HIq(2); R dl ). Since q G H 1 (P;IR dl ), it is sufficient to check 
q z G Ml(V;R dl ) for each z G dT> (recall (OSjl ). Since v z = q z o $ G Hj(R d ;M dl ), by 
virtue of Lemma 14.51 we get q z G R dl ). The proof is complete. □ 

Remark 4.1. The constant C appearing in estimate f !4.13|) does not depend on |o x | and 

Wx\- 

Proceeding identically to the proof of Proposition I4.6[ we can prove the following 

Proposition 4.7. Let the conditions of Theorem \ 3.1\ be satisfied. In addition, assume 
that the function pair (p,q) G T-L 2,1 (V) is a strong solution of BSPDE (11.11) and (II. 3p . 
Then (p,q) G $) 2 ' l (V), and there exists a constant C only depending on K,p ,n,T and 
the function 7 such that 

\\(p,q)\\h. H v)<C{\\\F\\\l v +Eml v ). (4.26) 
Now we use the standard method of continuation to prove Theorem 13.11 

Proof of Theorem \3.1[ The uniqueness of the strong solution to equation (II. ip is an im- 
mediate consequence of estimate (14.261) . Now we define 

Co = a ij (t, 0) Dij + b% x)Di - c(t, x), M k = cr ik (t, 0) A + u k (t, x), 
£1 = a ij (t,x)D ij + b l (t 1 x)D i - c(t,x), M\ = a ik (t 1 x)D i + u k (t,x). 

For A G [0, 1], define 

C x = (1 - A)£ + AA, M k x = (l-\)M k + \M k . 

Consider the following equation 

dp=-(C x p + M k q k + F)dt + q k dW k , P\ x&dv = 0, P\ 1=T = <P- (4-27) 

It is clear that the coefficients of equation (I4.27P satisfy the conditions of Theorem 13.11 
with the same K, k and 7. Hence a priori estimate (14.261) holds for equation (14.271) for 
each A G [0, 1] with the same constant C (i.e., independent of A). 

Assume that for some A = A G [0,1], equation (I4.27P is solvable, i.e., it has a 
unique strong solution (p, q) such that (p, q) G $) 2,1 (V) for any F G H°(T>) and any 
<p G L 2 (tt, J^, HHP)). For other A G [0, 1], we can rewrite (Q7II as 

dp = -{C Xo p + M\ Q q k + (A - A Q ) [(& - £ )P + (M k - M k )q k ] + F}dt + q k dW k . 
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Thus for any (u,v) G "% 2,1 (T>), the equation 

dp = -{C XoP + M k Xo q k + (A - A ) [(A - Cq)u + {M\ - M k )v k ] + F}dt + v k dW k , 

with the boundary conditions p\t=r = (ft and p\ x ^&v = 0, has a unique strong solution 
(p,q) such that (p, q) G i3 2,1 (T>). Then we define the operator 

A : H 2,1 (V) -> H 2,1 {V) 

as follows: 

A(u,v) = (p,g). 

Note that A(it, f) G -fj 2,1 (P). Then from estimate (14.261) . we can easily obtain that for 
any (u h Vi) G H 2 > l (V), i = 1,2, 

1 1 - U 2 ,«i - «z) ||wa,i(j3) < ||-A(Wl-W2,Ul-'W2)||Ja, 1 (jj) ( 4 - 2 8) 

< C|A- A |||(ui-m 2 ,«i -^Hwa,!^)- (4.29) 

Recall that the constant C in (I4.28[) is independent of A. Set 9 = (2C)~ l . Then the 
operator is contraction in 1-L 2,1 {T>) as long as |A — Ao| < 0, which implies that equation 
g2ZD is solvable if |A - A | < 9. 

Equation (14. 2 7ft is solvable for A = in view of Proposition 14.61 Starting from A = 0, 
we get to A = 1 in finite steps, and this finishes the proof of solvability of equation (11.11) . 
The proof of Theorem 13.11 is complete. □ 

5 Some applications 

5.1 A comparison theorem 

The comparison theorem plays an essential role in the theory of PDEs and BSDEs. Ma 
and Yong [15] gives some comparison theorems for strong solutions to the Cauchy problem 
of degenerate BSPDEs by Ito's formula, which are improved by Du and Meng [5] under 
the super-parabolicity condition. In this subsection, we prove the following comparison 
theorem for the strong solution to BSPDE (11.11) and (II. 3p in the general C 2 domain. 

Theorem 5.1. Let the conditions of Theorem \3.1\ be satisfied, and (p,q) be the unique 
strong solution to BSPDE (II. ip and (II. 3p . Suppose for any (cu,t), F(t, •) > and > 0. 
Then p(t, •) > a.s., for every t G [0,T]. 

The proof of Theorem 15.11 needs the following lemma. In what follows, we denote 
a~ = -(a A 0) for a G R. 

Lemma 5.2. Let the conditions of Theorem \3.1\ be satisfied. In addition, assume that a x 
and a x are bounded (by a constant L). Let (p,q) be the strong solution of equation (II .2p . 
Then for some constant C , 

E [ [p(t,x)-] 2 dx Kc^-'Ae [ [4>(x)-] 2 dx + E f T [ [F(s,x)-] 2 dxds\. (5.1) 
Jv I Jv Jt Jv ) 
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Proof. Define the function h : R — >■ [0, oo) as follows: 

{r 2 , r < —1, 

(6r 3 + 8r 4 + 3r 5 ) 2 , -1 < r < 0, 
0, r > 0. 

Then /i is C 2 and 

/i(0) = h'(0) = h"(0) = 0, h(-l) = 1, /i'(-l) = -2, = 2. 

For any e > 0, define /i e (r) = e 2 h(r/e). The function /t e has the following properties: 
\\mhJr) = (r~) 2 , lim h! (r) = — 2r~, uniformly; 

£->0 £-5-0 

|/i"(r)| < C, Ve>0,rGf; lim/i'Yr) = ( )?' r < °' 

e^o £ w [ 0, r > 0. 

Since a x and <t x exists and they are bounded, equation (II .ip can be written into the 
divergence form. Then applying Ito's formula for Hilbert-valued semi-martingales (see 
e.g. [H Page 105]) to h £ {p{t, ■)), and from Green's formula, we obtain that 



E \ h £ (cj)(x))dx — E / h £ {p{t,x))dx 
Jv Jv 

= Ej J ^-K( P )D t (a l W 3 p + a lk q k )-h' £ ( P )[(b l -D 3 a i:i )D lP 
-cp+ (y k - D i0 - ik )q k + F] + ^K(p)\q\ 2 ^dxdt 

= E it L {k h '^ P) ( 2alW *P D iP + Zc- lk q k D lP + \q\ 2 ) 

- h' e (p) [(b l - D ja ij )D iP -cp+ {v k - D i0 - lk )q k + F] ^dxdt. 
Setting e — > 0, by Lebesgue's Dominated Convergence Theorem, we have 

E i [<j)(x)~fdx - E i [{p{t,x)-} 2 dx 
Jv Jv 

= E I [ 1 {P < 0} • { (2a ij D iP D jP + 2a ik q k D iP + |g| 2 ) 
Jt Jv 

- 2p[(b i - D ja ij )D iP -c P + (u k - Dia ik )q k + F] }dxdt. 
For positive numbers S, Si, we have 

x 

2a ij D lP D jP + 2a lk q k D iP + \q\ 2 > 2a ij D iP D jP - (1 + 5)|a i Ap| 2 + 7^l<?| 2 

> { .2SK + (l + S)K}\D P \ 2 + -^-\q\ 2 , 

l + o 

- P [{U - D^)D iP — c P + (u k - D^ k )q k ] > -Si(\D P \ 2 + \q\ 2 ) - C{K, L)6?\pf 
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Taking 5 and 5i small enough (such that 5i = min{— 25K + (1 + S)k, j^} > 0), we have 

E I [<p(x)~} 2 dx — E I [(p(t,x)-} 2 dx 
Jv Jv 

>E f [ 1 {P < 0} - [-C(n,K,L)\p\ 2 -2pF]dxdt 

Jt Jv 

>E [ [ [-C(k,K,L)\ p -\ 2 -2p-F~]dxdt 

Jt Jv 

>E I [ [-C(n,K,L)\p-\ 2 -\F-\ 2 ]dxdt, 

Jt Jv 

and this along with the Gronwall inequality implies the desired inequality (15. ip . □ 

Proof of Theorem \5.1[ Due to Assumption 13.31 we can construct (e.g., by the standard 
technique of mollification) two sequences a n and a n with bounded first derivatives in x, 
which converge uniformly (w.r.t. (uj,t,x)) to a and a, respectively, with a n ,a n satisfying 
the same assumptions as a, a (with k 2 instead of k). Then, in view of Theorem 13. 1[ the 
following BSPDE (for each n) 



dp n = -(a%D ijPn + VDiPn - cp n + a l n k D iq * + v k q k n + F)dt + q k dW t 



k 

t -i 



PnlxedV = 0, Pn\t=T = 4> 

has a unique strong solution (p n , q n ) 6 I 2 x I 1 , such that 

llb»lll2 + lll?n|||?+£ sup \\ Pn (t,-)\\l<C( HI +EU\\ 2 ), (5.2) 

0<t<T 

where the constant C only depends on K, p , k, T and the function 7, and does not depend 
on n. It is easy to check that the function pair (p—p n , Q—Qn) satisfies the following BSPDE 



du = - (a ij D ijU + UDiU -cu + a ik D { v k + v k v k + F n )dt + v k dW t 



k 

t ) 



u\ xe av = 0, u\ t =T = 0, 
with u and v being the unknown functions, where 

F n = (a« - 4)D l3 Pn + (a 1 * - < k )D t q k . 
In view of (15. 2p and keeping in mind the uniform convergence of a n and o~ n , we have 

|||-Fi|||o "^0 as n —7- 00, 
and this along with estimate (I3.3P implies that for every t e [0, T], 

p n (t, •) p(t, ■), strongly in L 2 (Q x V). 

On the other hand, it follows from Lemma [5.21 that p n (t, •) > a.s. for every t £ [0, T]. 
Hence we get p(t, •) > a.s. for every t £ [0, T]. The proof is complete. □ 
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5.2 Semi-linear equations in C 2 domains 

Consider the following BSPDE: 

dp(t,x) = - [a ij (t,x)D ij p( y t,x) + a ik (t,x)Diq k (t,x) + F(t, x, p(t, x) , q(t, x))]dt 

+ q k (t, x)dW k , (t, x) G [0, T) x V; 
p(t,x) = 0, te [0,T], x G dV; 

p(T,x) = <f>(x), x G V. 

(5-3) 

Such a BSPDE is associated to a FBSDE in a similar way as shown in Tang [2T] . 

Assumption 5.1. Forany(w,v) G (H^(V)nH 2 (V)) xH\V; R dl ), the function F(t, x, u, v) 
is predictable as a function taking values in H°(T>). The function F is continuous 
in (u, v). Moreover, for any e > 0, there exists a constant such that, for any 
(ui,Vi) G (^oOP) n ^ 2 (^)) x H^V-tR* 1 ), i = 1,2 and any we have 

-,Mi,fi) - F(t, -,m 2 , ^2) ||o,x? < e{\\ui - u 2 \\ 2 ,v + \\vi - v 2 \\i,t>) 

+K £ (\\u 1 -u 2 \\ ,v+ IK -v 2 \\ ,v)- (5.4) 

Remark 5.1. Take F(t,x,u,v) = b l DiU — cu + is k v k , and equation ( 15. 3D becomes BSPDE 
(11.11) and (11.31) . Under the boundedness of b, c and u, we can easily check condition (15. 4p 
by the interpolation inequality. On the other hand, we shall see that the strong solution 
of BSPDE (El) belongs to Hl'^V) (c l-L 2 ' l {V)), which allows us to discuss equation (Oj) 
with the coefficients c, 1/ blowing up near the boundary of T>. However, we prefer not to 
purse this issue any further in this paper. 

Then we have the following 

Theorem 5.3. Let the functions a and a satisfy Assumptions \3.S\ and \3.3\ Let Assurap- 
tions \3.1\ and \5.1\ be satisfied. Suppose 

F(-, ■, 0, 0) G E°(V), G L 2 (Q, & T , Hl{V)). (5.5) 

Then BSPDE (15.31) has a unique strong solution (p,q) G -fj 2,1 (2}) such that 

||(p,g)||^ ) <C'(|||F(-,-,O,O)|||^+ J E||0||^), (5.6) 

where the constant C only depends on K,p , k,T, the functions 7 and K e . 

Proof. Step 1. First, we prove the a priori estimate (15. 6ft for the strong solution of equation 

(ED. 

From estimate (13.31) . there exists a constant C depending only K, po, k,T, and the 
function 7, such that 



2 



<C(\\\F(.,.,p,q) +EU\\l v ). 



On the other hand, in view of Assumption 15. H we have 

III F(; ;P, 9) < 2 III F(; ;P,q)~ ; 0, 0) +2 ||| F(-, ■, 0, 0, mv 
< 4e 2 ( HI p + HI q \\\l v ) + AK 2 ( HI p \\\l v + HI q |||^ ) + 2 ||| F(; •, 0, 0) |||^ . 
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Therefore, taking e = (8C) 2 , we have 

< <7( HI F(-, 0, 0) \\\ 2 >v +EU\\l v + ||b HIS.2, + III Q \\\l,v ) • (5-8) 

Proceeding similarly as in the proof of Proposition 14.61 we can remove the last two terms 
in the last inequality. Indeed, applying Ito's formula for Hilbert-valued semi-martingales 
to [pit, -)] 2 , we have 

\\p(0,-)\\lv=¥\\lv + 2 f T I p(a l W t3 p + a lk D iq k + F(t,x,p,q))dxdt 

Jo Jv 

\\q(t,-)\\l v dt-2 [ [ pq k dxdW t h . 
Jo Jv 

Taking expectations, using the Cauchy-Schwarz inequality, and keeping (15. 7p in mind, we 
have for any 5 > 

\h\\\lv < EU\\l v + 2E [ [ p(a ij D ijP + a ik D lQ k + F{t,x,p,q))dxdt 

Jo Jv 

< E\\4>\\l v + 5(\\\p \\\1 >V + HI q \\\l v + HI F(-, -,p, q) \\\l v ) + C(S, K) \\\ p\\\\ v 

< EU\\l v + 6( HI p HI 2 , P + (I) q \\\l v ) + C(5, K) ( (I) p \\\l v + (I) F(-, •, 0, 0) ). 

Taking 5 small enough and repeating (15.81) . we have 

\\(p, q)\\l^ {v) < C( HI F{; 0, 0) ||| 2 iP +EU\\l v + \\\p \\\l v ), 

where the constant C only depends on K, p , k, T, and the functions 7 and K £ . Using the 
Gronwall inequality, we obtain a priori estimate (15.61) . 

Furthermore, from a similar argument as above, we can prove the uniqueness of the 
strong solution of BSPDE Q . 

Step 2. We use the method of continuation to prove the solvability of BSPDE (15.31) . 
For each AG [0, 1], we consider the equation 

dp = - [a ij D ijP + a ik D iq + XF(t, x, p, q)] dt + q k dW k , p\ xe&D = 0, p\ t=T = 0. (5.9) 

It is clear that the function XF satisfies Assumption 15. II with the same K £ as F, and then 
equation (15.91) has a priori estimate (I5.6P for each A with the same constant C . 

Assume that for A = A e [0, 1], BSPDE (15.91) has a unique strong solution (p,q) G 
$) 2 ' 1 {V), for any <\> G L 2 (fi, & T , H£(V)) and any F G M°(V) satisfying Assumption O 
and condition (15.51) . For other A G [0, 1], we can rewrite (14.271) as 

dp = -{a ij D ijP + a ik D iq + X F(t, x,p, q) + (A - X )F(t, x,p, q)}dt + q k dW k . 

Thus for any (u,v) G T-L 2,1 (V), the equation 

dp = -{a ij D ij p + a ik Diq + X F(t,x,p,q) + (A - X )F(t, x, u, v)}dt + v k dW k , 

with the boundary conditions p\t=r = and p\ x€ ov = 0, has a unique strong solution 
(p,q) G S) 2,1 {V). Then define the operator 

A : n 2A (V) T-L 2,1 (V) 




24 



as A(u,v) = (p,q)- Note that A(u,v) G S) 2,1 (V). Proceeding similarly as in Step 1, we 
can easily obtain that for any (ui,Vi) G T-L 2,1 {T>), i = 1, 2, 

- u 2 , v 1 - v 2 ) \\n^(v) ^ \\ A ( u i ~ u 2, Vi - v 2 ) \\%2a {v) (5.10) 

< C|A - A | || (til - U 2 ,Vi - V 2 )\\ 2 H 2,i {v y (5.11) 

Recall that the constant C in (15.101) does not depend on A. Set 9 = (2C) -1 . Then the 
operator A is a contraction in 1-L 2,l {T>) as long as |A — Ao| < 0, which implies that (15. 9p is 
solvable if |A — A | < 6. 

The solvability of equation (15. 3ft for A = has been given by Theorem 13.11 Starting 
from A = 0, we can reach A = 1 in finite steps, and this finishes the proof of solvability of 
equation (15.31) . □ 
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